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ABSTRACT 



Context. The bulge dominated galaxy NGC 7814 provides one of the strongest dynamical tests possible for Modified Newtonian 
Dynamics (MOND). Spitzer 3.6 \im photometry fixes the bulge parameterisation and strongly constrains the properties of the sub- 
dominant stellar disk. Furthermore, the distance is known to better than 5%, virtually eliminating it as a free parameter. The rotation 
curve is easily measured, since the H I (and stellar) disks are edge on, and both the receding and approaching sides agree very well. 
Aims. To explore the agreement between the model and observed rotation curves in MOND given that the only two free parameters 
available are the mass-to-light ratios of the bulge and disk. 

Methods. We use a grid based MOND Poisson solver that accurately solves for the MOND gravity and produces our model rotation 
curves from a given mass distribution. The input to the Poisson solver is a 3D distribution of N particles which is generated from 
modelling the observed distribution of stars and gas in the galaxy. 

Results. By ensuring a superior fit to the radial surface brightness profile than previous works, by virtue of a double Sersic fit to the 
bulge, we were able to produce excellent fits to the rotation curve with typical values for both mass-to-light ratios. 
Conclusions. The model rotation curve of a mass distribution in MOND is extremely sensitive to the bulge-disk decomposition and 
even slight deviation from the observed mass distribution can produce large differences in the model rotation curve. 
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1. Introduction 

As H I studies of galaxies become increasingly more deep 
and resolved, both spatially and in terms of their velocity, the 
ability to constrain Modified Newtonian Dynamics (MOND; 
Milgrom 1983, but for a recent review see Famaey & McGaugh 
2012) grows. Arguably, MOND surpasses theories that require 
galactic dark matter in explaining the dynamics of galaxies 
mainly due to its greater predictiveness. 

The premier tool for measuring the dynamics of galaxies is 
the HI rotation curve. Once one is in possession of a galaxy's 
distance and luminous mass distribution, one should be able to 
precisely determine that galaxy's expected rotation curve and 
test the theory by comparing with the observed rotation curve. 
Unfortunately, uncertainties often get in the way of a straight- 
forward comparison between the measured and predicted rota- 
tion curve, whether be it in terms of the distance uncertainty, or 
more commonly from the so-called mass-to-light ratios of the 
stellar components. 

The mass-to-light ratio (M /L) of a stellar component is the 
unknown relationship between the luminosity, of that ensemble 
of stars in a certain waveband, and its mass. If the precise lu- 
minosity function of stars for the stellar component was known, 
then we could exactly determine the M/L, but not even the ini- 
tial mass function is known and multiple epochs of star forma- 
tion occurring at different locations within the galaxy may cloud 
the result from a simple single epoch. 

To combat this, state of the art photometry of galaxies is per- 
formed at infrared wavelengths such as the 2 fim-band (used by 
the 2MASS survey; Skrutskie et al. 2006) and 3.6 fim (used by 



the SINGS survey; Kennicutt et al. 2003). These wavelengths 
are more sensitive to the lower mass stars, which contribute the 
bulk of the mass, and are less sensitive to dust extinction than are 
optical wavelengths. Incorporating information about the colour 
of the stellar component enables us to date it and make our esti- 
mation of the M I L. 

As alluded to above, distances can also inject a large uncer- 
tainty, however, recent measurements based on locating the tip 
of the red giant branch in the galaxy's colour-magnitude diagram 
can facilitate a distance estimation with as little as five per cent 
error (e.g. Ferrarese et al. 2000). When one analyses the light 
distribution in a disk galaxy, if it is slightly inclined from being 
totally edge on, it is impossible to accurately determine the ver- 
tical distribution of stars and gas so one also has some freedom 
in fitting for the scale heights of those distributions. 

In addition to the uncertainties associated with the vertical 
mass distribution, the Mj Ls of the various stellar components 
and distance, there is also the inherent uncertainty associated 
with how MOND transitions between the strong and weak grav- 
ity regimes. Strong and weak are relative to the new, presum- 
ably fundamental, constant of acceleration, a , that encapsulates 
MOND. This also has uncertainty and is taken to be around 
a a w 3.6(kms _1 ) 2 pc _1 . This, preferably smooth, transition 
can take on a wide variety of forms and can substantially alter the 
predicted rotation curve at intermediate values of acceleration. 

Fraternali et al. (2011, hereafter FSK11) have recently pre- 
sented a comprehensive analysis of a bulge dominated galaxy - 
NGC 7814 - with an edge on disk that has a symmetrical rota- 
tion curve (receding and approaching) which can provide an ex- 
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cellent test of MOND, perhaps one of the strongest tests to date. 
For one, the distance is measured to better than 5 per cent using 
the tip of the red giant branch, the galaxy has 3.6 fim photometry 
and since it is edge on, it has an accurately determined vertical 
stellar distribution. The only tunable parameters that can make 
a significant difference to the fitted rotation curve are the M/Ls 
of the stellar components. The reason we talk about stellar com- 
ponents is because NGC 7814 is a bulge dominated galaxy, but 
it has a stellar disk in addition to an H I disk. 

NGC 7814 is one half of a pair of galaxies, the other being 
NGC 891, originally considered by van der Kruit (1983, 1995). 
These two are akin to the pair NGC 2403 and UGC 128 later 
considered by de Blok & McGaugh (1998) in the sense that it is 
of fundamental importance whether galaxies with the same total 
mass, but different mass distribution have appropriate rotation 
curves in MOND. 

Angus et al. (2012) presented a Poisson solver for Milgrom's 
QUasi-Linear formulation of MOND (QUMOND; Milgrom 
2010). It takes the input from any 3D distribution of particles 
and accurately calculates the QUMOND gravity everywhere. 
This is important because strictly the typical algebraic relation 
of MOND used by the majority of researchers is only valid in 
spherical symmetry. The code was shown to be effective at fit- 
ting the rotation curves of a small sample of the THINGS survey 
galaxies (de Blok et al. 2008). 

Here we seek to discover whether a detailed analysis of the 
available data on NGC 7814 and using our sophisticated tool 
to generate the MOND gravity from the source 3D density dis- 
tribution can assist an excellent fit to the rotation curve with 
physically motivated parameters i.e. reasonable M/Ls for both 
the bulge and the disk. In this paper the acceleration constant 
of MOND is always assumed to be a = 3.6( kms _1 ) 2 pc _1 
and the interpolating function is assumed to be v(y) = 0.5 + 
0.5-^/ 1 + 4/y, where y = v$iy l and <j>jv is the Newtonian po- 
tential. 

2. Bulge-Disk Decomposition 

Fits to an observed surface brightness with even mild central 
deviations - especially in the case of a bulge - can have a severe 
impact on the potential agreement between the rotation curve 
and MOND. Therefore, we made our own careful bulge-disk de- 
composition for NGC 7814. The primary piece of data we used 
was the radial surface brightness along the major axis. We took 
the data from Fig 6, right hand panel, of FSK11 and used the 
formula 2.5logioI(R) = 21.572 + 3.24 + ^3. 6pm to convert the 
surface brightness in magnitudes per squared arcsecond to pro- 
jected luminosity £o,3.6 M mPc~ 2 > where the 3.24 corresponds to 
the absolute magnitude of the Sun at 3.6/xm. 

Using a single Sersic profile it is not possible to obtain a good 
match to the surface brightness profile both for R < 0.5 kpc and 
0.5 < R < 1 kpc, where the bulge is the dominant compo- 
nent. Additionally, the stellar disk is the dominant component of 
surface brightness in the disk plane for R > 1 kpc (see Fig Q]), 
but it is also not possible to fit the surface brightness profile for 
1 < R < 4.5 kpc and R > 4.5 kpc with a single exponential 
disk. 

Therefore we fitted a double exponential model for the disk 
and double Sersic density profile for the bulge i.e. one described 
by a certain Sersic profile (or exponential disk) up to a given 
radius and then described by another beyond that radius. The 
combined, bulge and disk respectively, projected luminosity 
I(R, z) = I b {R, z) + I d {R, z) was thus 



I{R,z = Q) = I eea ;p(-7.67((i?/i? e ) 1/ns - 1)) 
+E d exp(-R/R d ) 

I(R = 0,z) = / e exp(-7.67((|z|/z e ) 1 /™= - 1)) 

+^ d sech 2 {\z\/z d ) (1) 



One combination of Sersic parameters (I ei , R ei , n Si ) is used 
only to fit the inner 0.5 kpc and another combination (I eo , R& , 
n So ) is used only beyond 0.5 kpc. Also one combination of disk 
parameters (S di , R di ) is used to fit the inner 4.5 kpc and another 
is used only beyond 4.5 kpc (T, do , R do ). We used the subscripts 
"i" and "o" refer to the inner and outer parameters, as per ta- 
ble Q] These parameters are not to be interpreted physically in 
terms of scaling relations. One can see that the former and lat- 
ter combinations of parameters bracket those of FSK11, in the 
sense that I ei is lower than the single value of the Sersic fit per- 
formed by FSK1 1 and I Eo is larger. All fit parameters are given 
in table 1 . 

One can see in Fig Q] that the radial surface density of the 
galaxy is better represented in our model (solid green line) than 
the single Sersic, single exponential disk profile fit by FSK1 1 
(dashed red line), especially below 0.3 kpc (right hand panel) 
and between 2 and 4.5 kpc (left hand panel). Beyond R = 
4.5 kpc, our exponential disk scale-length changes to a lower 
value and the surface brightness of our model goes through the 
centre of the data points, as does the FSK1 1 model. The dot- 
ted blue line of Fig Q] shows the surface brightness of our model 
with only a bulge. This is simply to emphasise that the surface 
brightness for R > 1 kpc is mainly contributed by the disk, not 
the bulge, even though the bulge still contributes the bulk of the 
gravity at all radii. 

In order to constrain the stellar disk scale-height and the axis- 
ratio of the bulge, q = z e /R e (both of which are constant with 
radius), we attempted to match the contours of surface bright- 
ness taken from FSK1 1 Fig 6 top panel as well as the z-direction 
surface brightness profile at R = kpc. The 2D contours of 
surface brightness are equally well represented in both models 
(see Fig 01, which shows that both sets of disk parameters are 
acceptable. In fact, many correlated combinations of central sur- 
face density and scale-length will allow agreement, however, the 
scale-height cannot differ greatly from our 0.42 kpc (or FSK 1 1 's 
0.44 kpc) or the outer contours will cease to agree. 

So that we can solve for the rotation curve of our model 
galaxy, we must produce a 3D model of the galaxy. In order 
to minimise the complexity, we only allow for the z-direction 
variation in terms of the axial ratio, q, rather than a fitting for a 
whole new set of Sersic parameters. Nevertheless, the data are 
relatively well adhered to (see Fig [2). 



2.1. Total luminosity 

Given this special case of a completely edge on disk, the lu- 
minosity must be inferred by integrating the projected luminos- 
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ity 

/>4.5 kpc 

L d = [/ 2Z dz exp(-R/R dt )dR + 
Jo 

/•OO 

/ 2Z do exp(-R/R do )dR] x 

J 4.5 kpc 

/oo 
sech 2 (z/z d )dz] 
-oo 

= 8 x 10 9 L a3 . 6A1 m- (2) 

The bulge luminosity must be integrated numerically and is 
found to be — lO 11 L0 3.g Aim , which is 43% higher than 
FSK11 bulge luminosity. 

It is worth remarking that these parameters, both in the text 
above and table Q] are all projected parameters. R is the major 
axis which is in the horizontal plane of the sky and z is the minor 
axis and is in the vertical plane of the sky. 

2.2. Distance and Rotation Curve 

The distance to NGC 7814 used for presenting the analy- 
sis in FSK11 was 14.6 Mpc, which is well within the errors of 
the distance derived by Radburn-Smith et al. (2011), as part of 
the GHOSTS survey, which is 14.4t° g Mpc (cf. their table 5). 
For ease of comparison, the distance is always assumed to be 
14.6 Mpc in this paper, it is not a free parameter. 

The rotation curve we use for comparison with the models is 
the mean of the approaching and receding rotation curves in the 
position velocity diagram presented in Fig 3 of FSK1 1 (private 
communication with the authors of FSK1 1). 

3. Our model 

We generated a 512 3 particle realisation of NGC 78 14 where 
we dedicated half of the particles to the bulge and a quarter each 
to the H I and stellar disks. The luminosity profile was strictly 
adhered to using the surface brightness models advocated above, 
as was the H I surface density, found in FSK1 1 Fig 7 right hand 
panel, by making an Abel transformation from our analytical 2D 
models to numerical 3D models, which we stored. We assume 
the galaxy is symmetrical around the minor axis. 

For distributing the particles of the bulge, we used a standard 
rejection technique: we drew random 3D radii from an exponen- 
tial distribution and simply interpolated through the data points 
of our numerical 3D model (created by the Abel inversion) to 
find the 3D density at that given radius. Then we decided to keep 
the particle if the next random number was less than the ratio of 
the value of the density at that radius of the numerical density 
and the analytical one. We then randomly assigned x, y and z 
spherical coordinates but multiplied the spherical z coordinate 
by the axis ratio q to create the flattening. For the disks we dis- 
tributed the particles in z according to a sech 2 (z) profile: the 
gas with a 0. 1 kpc scale-height and stars using the scale-height 
given in table Q] The x and y coordinates were chosen as for the 
bulge, except our analytical distribution was 2D; meaning the 
cumulative distribution function differs from the 3D case. 

Using the QUMOND Poisson solver code presented in 
Angus et al. (2012), we solved for the gravity in the midplane 
of the disk and used this to produce the model values of the ro- 
tation speed as a function of projected radius, V mode i(R). 

Our task then was to explore the parameter space of 
the free parameters in order to elucidate the agreement of 



MOND with the observed rotation curve. The only param- 
eters free to be varied were the mass-to-light ratios of the 
bulge and disk. We parameterise the goodness of fit as x 2 = 

( ^g^^ ) 2 , where n= 18. 
4. Results 

In Fig [4] we show contours of goodness of fit, X 2 /n for the 
2D parameter space defined by the mass-to-light ratios of the 
bulge and disk (x and y axes respectively). We find that the best 
fits (x 2 /n < 0.35) occur for values from M/Ld = 0.6 to 1.1 
and that values as low as 0.3 (and high as 1.4) can still produce 
excellent fits with x 2 /n < 0.45. To achieve x 2 / n < 0-45 the 
mass-to-light ratio of the bulge, M/Lb, is required to be in a 
significantly tighter range 0.58 to 0.685. 

In Fig [5] we plot the rotation curves of four models against 
the data points measured by FSK1 1 . The black lines use our pa- 
rameterisation of the bulge and disk and use mass-to-light ratios 
of (0.66,0.7) and (0.685,0.3) for the solid and dashed lines re- 
spectively. 

4. 1. Comparison with other work 

FSK11 performed a bulge-disk decomposition which sepa- 
rated the stellar distribution into a bulge and thin disk. To test 
MOND, they used the velocity profiles derived for each com- 
ponent (i.e. they solved Poisson's equation to convert surface 
brightness to gravity and hence circular speed) and they found 
that agreement between the model rotation curve and the mea- 
sured one could only be found if M/ Ld = 4.6; an unfeasibly 
large disk mass-to-light ratio. In Fig[TJthe red dashed line corre- 
sponds to the FSK1 1 fit to the surface brightness, which is not as 
good a fit as our model to the inner 4.5 kpc. This is simply be- 
cause NGC 7814 has a relatively complicated surface brightness 
profile and a double Sersic bulge and double exponential disk is 
required to trace the profile over the full range. 

In order to compare the surface brightness models of FSK1 1 
we need to know their parameterisations. The difference be- 
tween projected and deprojected parameters is never clearly dis- 
tinguished in FSK11, but we assume the "r e " they use is the 
projected Sersic radius. On the other hand, the surface density 
they give for the disk is almost certainly projected in z and de- 
projected along the line of sight (i.e. as seen from above the disk 
- face on, whereas our true perspective is that of an edge on disk), 
even though the disk scale-length given is presumably the pro- 
jected one. Therefore, using the luminosity of their disk, we re- 
derived the parameters for their model and have added them to 
table[T] 

Using these parameters we were able to follow the same pro- 
cedure as describe in the beginning of Sj3]to generate a 3D model 
of NGC 7814. We then fully explored the parameter space de- 
fined by the mass-to-light ratios of the two stellar components. 
The plain contours in Fig|4]are using FSK1 l's parameterisation, 
for which the best fit clearly requires a large Mj Ld of around 
2.65 and it must be larger than 1.5 in order to have a x 2 /n < 0.8. 
This is incompatible with the expected mass-to-light ratios of 
star forming regions in the 3.6 /im band, which should have val- 
ues lower than that of the bulge, which itself should have a value 
around 0.7. The best fit mass-to-light ratio for their bulge is 0.61- 
0.73, which is not too dissimilar to our own. 

In Fig [5] we have plotted the best fit rotation curve using 
the FSK11 model (red solid line), with M/L B = 0.66 and 
M/Ld = 2.65, against the measured rotation curve as well as 
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Fig. 2. Here we plot the projected luminosity against projected 
radius along the minor axis of NGC 7814. The data points come 
from FSK11 and their fitted model is the red dashed line. The 
green solid line is our model, which is not a fit, but rather is 
the projected luminosity profile fit of the major axis squashed 
by the axis ratio, q, which is determined from fitting the surface 
brightness contours of Fig [3] 
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Fig. 3. Here we show the fits to the surface brightness contours 
of NGC 7814 presented in FSK11. The models are the closed 
curves and the filled circles are the observed contours. The black 
lines correspond to our model and the red lines to that of FSK1 1 . 
The parameters of the two models are given in table Q] 



a fit using a more regular value for the disk mass-to-light ratio, 
M/Lb = 0.75 and M/L D = 0.7 (red dashed line). The best 
fit is, as expected, as good a fit as the solid black line (our best 
fit model), but the red dashed line shows that using a regular 
M/Le> leads to a very poor fit. 

We emphasise that it is the difference between our brightness 
profile for the bulge and that of FSK1 1 that allows us to reduce 
the mass-to-light ratio of the subdominant stellar disk, not the 
different parameterisation of our disk and theirs. 



5. Conclusion 

The bulge dominated galaxy NGC 7814 has an edge on stel- 
lar disk and a regular H I disk. The bulge-disk decomposition 
performed by FSK11 suggested that an impossibly large mass- 
to-light ratio for the stellar disk was required in order to fit the 
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Fig. 4. Here we show the contours of goodness of fit, x 2 /n for 
the 2D parameter space defined by the mass-to-light ratios of the 
bulge and disk (x and y axes respectively) found by fitting the 
rotation curve of NGC 7814. The filled contours are found using 
our parameterisation of the bulge and disk and the plain contours 
are using FSKll's parameterisation. Both sets of contours use 
the same levels and both use the simple v function of QUMOND. 
The central contour has a value x 2 jn = 0.35 and the subsequent 
contours are 0.45, 0.8, 2 and 4. 
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Fig. 5. Here we show the rotation curves of four fits to the mea- 
sured rotation curve of NGC 7814 (data points from FSK11). 
The black lines use our parameterisation of the galaxy and the 
red lines use that of FSK1 1. The bulge and disk mass-to-light ra- 
tios are (0.66,0.7), (0.685,0.3), (0.66,2.7) and (0.75,0.7) for the 
linetypes black solid, black dashed, red solid and red dashed re- 
spectively. 



rotation curve of the galaxy. Here we used two Sersic profiles to 
model the projected luminosity of the bulge such that one pro- 
file is used below 0.5 kpc and the other is used beyond 0.5 kpc. 
This allows a superior fit to the observed surface brightness, es- 
pecially below 0.3 kpc and between 2 and 4.5 kpc. We then de- 
projected the projected luminosity models for bulge and disk, 
modelled the galaxy with a 3D particle distribution and used the 
code of Angus et al. (2012) to compare the observed and mod- 
elled rotation curves across the 2D parameter space comprising 
the mass-to-light ratios of the bulge and disk. We made the same 
analysis using the models of FSK1 1. 
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range \/R 2 + z 2 /q 2 < 0.5 kpc 


y / R 2 +z 2 /q 2 > 0.5 kpc 




R < 4.5 kpc 


R > 4.5 kpc 








model 


I ei -Re; n Si 




q 


Srf; Rdi 


Sd Rd a 


Zd 


Lb Ld 




Lqpc~ a kpc 


Lqpc~ 2 kpc 




Lqpc~ 2 kpc 


Lqpc~ 2 kpc 


kpc 


10 iu L Q 10 1U L Q 


A12 


60.05 3.82 1.747 


9000 0.578 5.747 


0.61 


900 6.0 


2250 2.7 


0.42 


10.0 0.8 


FSK11 


1120 2.16 4 


1120 2.16 4 


0.61 


1200 4.26 


1200 4.26 


0.44 


7.0 0.84 



Table 1. Here we list the parameters of the bulge and disk in our fitted model (referred to as A12) and the FSK1 1 fit. Note the ranges 
of applicability for the A12 model. The subscripts "i" and "o" refer to the inner and outer fitted parameters, both for the bulge and 
disk. All luminosities are in the 3.6 [im band. 




Fig. 1. In the left hand panel we plot surface brightness against projected radius along the major axis of NGC 7814. The data points 
come from FSK1 1 and their fitted model is the red dashed line. Our model, using the combined luminosity of the bulge and disk, 
is the green solid line and the blue dotted line is the surface brightness of our fitted bulge only. The right hand panel shows the 
projected luminosity against projected radius along the major axis. The fitted models are the same as in the left hand panel. 



It was found that with the FSK1 1 model, excellent fits to the 
rotation curve could only be found using very large mass-to-light 
ratios for the disk (larger than 1 .5 in order to have a goodnes s of 
fit x 2 / n < 0-8). With our model it was found that x 2 /n < 0.45 
was achievable for disk mass-to-light ratios in the range 0.3-1.4 
and that the goodness of fit was most strongly constrained by the 
mass-to-light ratio of the bulge - which should be between 0.58- 
0.685 for the same level of goodness of fit. This mass-to-light 
ratio for the ancient population of stars in the bulge is precisely 
in the range expected by stellar population models (see Bell & 
de Jong 2001). 

Here we have re-analysed the required properties of one of 
the most pristine test galaxies of Modified Newtonian Dynamics. 
There is virtually no room for manoeuvring with respect to the 
galaxy's inclination, distance, structure and we keep all param- 
eters associated with MOND at the standard values used in the 
literature. We found that by forcing a more stringent adherence 
to the luminosity profile of the bulge than previously used, it is 
possible to, firstly, make an excellent fit to the rotation curve of 
the galaxy and, in doing so, isolate the two free parameters re- 
quired by the fit. These two free parameters, namely the mass-to- 
light ratios of the bulge and disk, are perfectly compatible with 
the predicted vales. 

Our results show that a proper comparison between observed 
rotation curves and the MOND expectations requires extreme 
care when performing bulge-disk decompositions, especially for 
the central regions. 
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